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ABSTRACT 

In this paper, we propose low-complexity robust adaptive beamform¬ 
ing (RAB) techniques that based on shrinkage methods. The only 
prior knowledge required by the proposed algorithms are the angular 
sector in which the actual steering vector is located and the antenna 
array geometry. We firstly present a Low-Complexity Shrinkage- 
Based Mismatch Estimation (LOCSME) algorithm to estimate the 
desired signal steering vector mismatch, in which the interference- 
plus-noise covariance (INC) matrix is estimated with Oracle Ap¬ 
proximating Shrinkage (OAS) method and the weights are computed 
with matrix inversions. We then develop low-cost stochastic gra¬ 
dient (SG) recursions to estimate the INC matrix and update the 
beamforming weights, resulting in the proposed LOCSME-SG algo¬ 
rithm. Simulation results show that both LOCSME and LOCSME- 
SG achieve very good output signal-to-interference-plus-noise ratio 
(SINR) compared to previously reported adaptive RAB algorithms. 

Index Terms — Covariance matrix shrinkage method, robust 
adaptive beamforming, low complexity methods. 

1. INTRODUCTION 

Several important applications of adaptive beamforming like wire¬ 
less communications, radar and sonar, microphone array processing 
have been intensively studied in the past years. However, under cer¬ 
tain circumstances, adaptive beamformers may suffer performance 
degradation due to short data records, the presence of the desired sig¬ 
nal in the training data, or imprecise knowledge of the desired signal 
steering vector. In order to address these problems, robust adaptive 
beamforming (RAB) techniques have been developed in recent years 
m El From a design principle point of view, the generalized side- 
lobe canceller, worst-case optimization (3, diagonal loading nm, 
eigenspace projection and steering vector estimation with presumed 
prior knowledge (Hill have been investigated. However, RAB de¬ 
signs based on these principles have some drawbacks such as their 
ad hoc nature, high probability of subspace swap at low SNR and 
high computational cost m 

Recent works have focused on design approaches that combine 
different principles together to improve RAB performance. Methods 
which jointly estimate the mismatched steering vector using Sequen¬ 
tial Quadratic Program (SQP) (8) and the interference-plus-noise co- 
variance (INC) matrix using a shrinkage method cni have been re¬ 
ported. Another similar approach which jointly estimates the steer¬ 
ing vector using an SQP and the INC matrix using a covariance re¬ 
construction method El, presents outstanding performance com¬ 
pared to other RAB techniques. However, their main disadvantage 
is the high computational cost associated with the optimization algo¬ 
rithms IIOIIIII and the matrix reconstruction process El,ll6|-[?]. 
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This paper proposes adaptive RAB algorithms with low com¬ 
plexity, which require very little in terms of prior information and 
shows a better performance than previously reported algorithms. 
Firstly, the steering vector of the desired signal is estimated using 
a Low-Complexity Shrinkage-Based Mismatch Estimation (LOC¬ 
SME) algorithm. An extension of the Oracle Approximating Shrink¬ 
age (OAS) method El is employed to perform shrinkage estimation 
of the cross-correlation vector between the sensor array received data 
and the beamformer output. The mismatched steering vector is then 
efficiently estimated without any costly optimization procedure in 
a low-complexity sense. Secondly, we estimate the desired signal 
power using the desired signal steering vector and the input data. 
We also develop a stochastic gradient (SG) version of LOCSME, 
denoted LOCSME-SG, which does not require matrix inversions or 
costly recursions. In particular, in LOCSME-SG the INC matrix 
from the input data is estimated using a Knowledge-Aided (KA) 
shrinkage 03 approach along with the computation of the beam¬ 
forming weights based on the estimated steering vector through SG 
recursions. The proposed LOCSME and LOCSME-SG algorithms 
circumvent the use of direction finding techniques for the interferers 
when obtaining the INC matrix and only require the angular sector 
in which the desired signal steering vector lies as prior knowledge. 

This paper is structured as follows. The system model and prob¬ 
lem statement are described in Section 2. The proposed LOCSME 
and LOCSME-SG algorithms are introduced in Sections 3 and 4, 
respectively. Section 5 presents and discusses the simulation results. 
Section 6 gives the conclusion. 


2. SYSTEM MODEL AND PROBLEM STATEMENT 


Consider a linear antenna array of M sensors and K narrowband 
signals. The data received at the ith snapshot can be modeled as 


x(() = A(0)s(i)-f n(i), (1) 

where s(i) £ are uncorrelated source signals, 6 = [0i, • • • , £ 

is a vector containing the directions of arrival (DoAs), A(0) = 

[a(0i) -b e, • • • , a(0K)] € matrix which contains the 

steering vector for each DoA and e is the steering vector mismatch 
of the desired signal, n(i) G assumed to be complex 

Gaussian noise with zero mean and variance cr^. The beamformer 
output is given by 

i/(i) = w^x(i), (2) 

where w = [wi, ■ • ■ , wm]"^ £ is the beamformer weight 

vector, where {■)^ denotes the Hermitian transpose. The optimum 
beamformer is computed by maximizing the signal-to-interference- 
plus-noise ratio (SINR) given by 


SINR = 


a1\w^ a.f 


w^Ri+„w' 


(3) 




where af is the desired signal power, Ri+n is the INC matrix and 
assume the steering vector a is known precisely (a = a(0i)), then 
problem (3) can be transformed into an optimization problem as 


H T-> 

minimize w Ki+nW 

w 

subject to w^a = 1, 


(4) 


which is known as the MVDR beamformer or Capon beamformer 

m . The optimum weight vector is given by MVopt = —. 

Since Ri+n is usually unknown in practice, it can be estimated by 
the sample covariance matrix (SCM) of the received data as 


R(*) = - ^x(fc)x"(fc), (5) 

i ^' 
k^l 

which results in the Sample Matrix Inversion (SMI) beamformer 
wsM/ = However, the SMI beamformer requires a large 

number of snapshots to converge and is sensitive to steering vector 
mismatches HOIII II . The problem we are interested in solving is 
how to design low-complexity robust beamforming algorithms that 
can preserve the SINK performance in the presence of uncertainties 
in the steering vector of a desired signal. 


3. PROPOSED LOCSME ALGORITHM 


where v = tr(S)/M and S = diag(xi/*). By shrinking S towards 
F CD and subsequently using it in a vector shrinkage form, taking 
into account the snapshot index, the result gives 

d(i) = p(i)diag(F(i)) + (1 - p(i))diag(S(i)), (9) 


which is parameterized by the shrinkage coefficient p{i). If we 
define D = diag(d), then the goal is to find the optimal value 

of p{i + 1) that minimizes the mean square error (MSE) of 

- - 2 

F[||D( i -b 1) — F(i) II ] in the ith snapshot, which leads to 

^ (1 - i|)4r(DWS*(t)) +tr(D(i))tr(D*(i)) 

^ (i -b 1 - ^)tr(D(i)S*(i)) + (1 - jg)tr(D(i))tr(D*(i)) ’ 

(10) 

where the derivation is shown in the Appendix and S(i) = diag(l(i)), 

^ i 

where l(i) = i Ti{k)y* (k), is the sample correlation vector 
k — 1 

(SCV). Alternatively equation (10) can be re-expressed in vector 
multiplication form and leads to the following 


^ (l-l|)d^(i)i(i) + tr(D(i))tr(D*(0) 

(t + l-^)d«(i)i(i) + (l-^)tr(D(t))tr(D*(i))’ 

( 11 ) 

As long as the initial value of p(0) is between 0 and 1, the iterative 
process in (9) and (11) is guaranteed to converge OH. Once the 
correlation vector d is obtained, the steering vector is estimated by 


In this section, the proposed LOCSME algorithm is introduced. The 
basic idea of LOCSME is to obtain a precise estimate of the desired 
signal steering vector and afterwards use it to estimate the desired 
signal power and to derive the recursion for the weight vector. The 
estimation of the steering vector is described as the projection onto 
a predefined subspace matrix of an iteratively shrinkage-estimated 
cross-correlation vector between the beamformer output and the ar¬ 
ray observation data. To obtain the INC matrix, in LOCSME we use 
OAS method to shrink the SCM in order to estimate the INC matrix. 

3.1. Steering Vector Estimation 


ai(i) 


Pd(i) 

l|Pd(i)ll2' 


( 12 ) 


3.2. Desired Signal Power Estimation 

This subsection will exploit a novel method to estimate the desired 
signal power a^. This can be accomplished by directly using the 
desired signal steering vector. Let us rewrite the received data as 

K 

x(z) = ai( 2 )si + ^ afcSfc + n(z). (13) 

k=2 


The cross-correlation between the array observation data and the 
beamformer output can be expressed as d = F{xj/*}. With as¬ 
sumptions that |a,rivv| <C |aiw| for m = 2, • • • ,K and that signal 
sources and that the system noise have zero mean while the desired 
signal is independent from the interferers and the noise, d can be 
rewritten as d = i?{cri^ai^wai -b nn^w}. By projecting d onto a 
predefined subspace which collects all possible information from 
the desired signal, the unwanted part of d can be eliminated. The 
prior knowledge amounts to providing an angular sector in which 
the desired signal is located, say \9i — 9^, 9i -b 9^]. The subspace 
projection matrix P is given by 

P = [ci,C2,--- ,Cp][ci,C2,--- ,Cp]", ( 6 ) 

where ci , • • ■ , Cp are the p principal eigenvectors vectors of the ma¬ 
trix C, which is defined by m 

C= J Si{9)a^{9)d9. (7) 

9i —9e 

We then employ the OAS shrinkage technique in order to achieve a 
more accurate estimation of d, so that it can help us to obtain a better 
estimate of the steering vector. Let us define 

( 8 ) 


Pre-multiplying the above equation by (i) and assuming ai (i) is 
uncorrelated with the interferers, we obtain 

af (i)x(i) = af (i)ai(i)si -b af (i)n(i). (14) 

Taking the expectation |af^(i)x(i)|^, we obtain 


|af (i)x(i)|' = E[(af (i)ai(i)si + af (i)n(i))* 

(af(i)ai(i)si-baf(i)n(i))]. (15) 

If the noise is statistically independent from the desired signal, then 
we have 


|af(i)x(i)f = |af(t)ai(t)nsif+ af(t)F[n(i)n"(i)lai(i), 

(16) 

where i5[n(i)n^(i)] represents the noise covariance matrix R„ 
which can be replaced by where is assumed known here 

for convenience, otherwise it can be easily estimated by a specific 
estimation method. Replacing the desired signal power |sif by its 
estimate (i), the desired signal power estimate is computed as 




|af(i)x(t)f-af(t)ai(t)a^ 

|af(t)ai(i)P 


(17) 


Equation (17) has a low complexity (0(M)) and can be directly 
implemented if the desired signal steering vector is well estimated 
and the noise level is known. 


F = vl, 








3.3. Estimation of the INC matrix 

In this subsection we describe a method to estimate the INC ma¬ 
trix that is based on the OAS method fT2l and used in LOCSME. In 
the OAS estimation of LOCSME, we need the SCM in (5) as a pre¬ 
liminary estimate for the INC matrix. Then we define Fq = vq!, 
where uo = tr(R)/M. By minimizing the MSE described by 
i?[||R(i) — Fo(i — 1) II I 1121 . the following recursion is employed: 

R(0 = po(i)Fo(i) -b (1 - Po(i))R(i), (18) 

- ^ (l-l|)tr(R(i)R(i))+tr^(R(i)) 

(* + l-^)tr(R(i)R(i)) + (l-^)tr2(R(*))’ 

(19) 

where /3o(0) must be initialized between 0 and 1 to guarantee con¬ 
vergence QIJ. To exclude the information of the desired signal from 
the covariance matrix of the array observation data, a simple sub¬ 
traction is considered: 

Ri+n(i) = R(i) - (3-i(i)ai(i)af (i). (20) 


Now, we employ the idea of KA shrinkage method ca to help with 
our INC estimation. By applying a linear shrinkage model for the 
INC matrix, we have 

Ri+„(i) = T]{i)Ro + (1 - ri{i))Ri+n(i), (24) 

where Rq is an initial guess for the INC matrix, ri{i) is the shrink¬ 
age parameter and rj{i) £ (0,1). Here the shrinkage parameter is 
expected to be adaptively estimated. Employing an idea of adaptive 
filtering [BJ. it is possible to set the overall filter output yf{i) equal 
to [R.i+„(i) ai(i)]^x(i) which is the linear combination of the out¬ 
puts from two filter elements which are yof{i) = [Roai(i)]^x(i) 
and yf{i) = [Ri+n(*)ai(i)]^x(i), which leads to 

yf{i) = v{i)yofii) + (1 - n{i))yf{i)- (25) 

To restrict the value of T]{i) equal to either 0 nor 1, a sigmoidal func¬ 
tion is employed: 

nii) = sgm[e(i)] = (26) 


3.4. Computation of Beamforming Weights 


where t(i) is updated as 


Eor the proposed LOCSME algorithm the beamforming weights are 
computed directly by 


Ri+^(i)ai(i) 

af(i)R-+i„(i)ai(i)’ 


( 21 ) 


which has a computational costly matrix inversion RY_^„(*). to 
reproduce the proposed LOCSME algorithm, whose complexity 
is equations (9),(11),(12) and (17)-(21) are required. In 

comparison to previously reported RAB algorithms in EllllQolIII] 
with costly online optimization procedures and complexity 0{M^) 
or higher, LOCSME requires a similar or lower cost. 


4. PROPOSED LOCSME-SG ALGORITHM 

In this section, the proposed LOCSME-SG algorithm is detailed. 
The aim is to devise a low-complexity alternative to LOCSME that 
is suitable for time-varying scenarios and implementation purposes. 
LOCSME-SG employs identical recursions to LOCSME to estimate 
the steering vector and the desired signal power, whereas the esti¬ 
mation of the INC matrix and the beamforming weights is different. 
In particular, LOCSME-SG employs a Modified Array Observation 
(MAO) vector to compute a preliminary estimate of the INC matrix, 
followed by a refined estimate with a low-cost shrinkage method. 

4.1. Estimation of the INC matrix 

In this subsection we present an extension of the KA shrinkage 
method (B) to estimate the INC matrix, which has much lower 
complexity than the one used in LOCSME. In LOCSME-SG, with 
the estimate of the desired signal power we subtract unwanted in¬ 
formation of the interferences out from the array received data in a 
low complexity vector form to obtain the MAO vector. Consider a 
simple substraction step as 

Xi+„(i) = x(i) - (Ti(i)ai(i). (22) 

Then the INC matrix can be estimated by 

Ri+n(*) = Xi+„(i)x" „(i). (23) 


e{i -b 1) = e{i) - - -|yo/(() - f//(()|^ 

(cTi +'?(*)) 

+ T^iiVofi'Q - W(*))y/(*)})»7(*)(1 - '?(*)). (27) 

where pe is the step size while is a small positive constant, q{i) 
is updated as 

q{i -b 1) = Aq(i)(l - Xq)\yof{i) - yf{i)\^, (28) 

where Xg is a forgetting factor (B) 


4.2. Computation of Beamforming Weights 

For the proposed LOCSME-SG algorithm, we resort to an SG adap¬ 
tive algorithm to reduce the complexity required by the matrix inver¬ 
sion. The optimization problem (4) can be re-expressed as 


minimize w^(i)(x(i)x^(i) 

w(i) 

subject to w^(i)ai(i) = 1. 


d-i(i)ai (i)af(i))w(i) 

(29) 


In order to compute the beamforming weights, we employ an SG 
recursion as given by 

dL 

w(i + l) = w(I)-M^, (30) 

whereE = w^(i)(x(i)x‘^(i)—(Ti(*)ai(i)ai^ (i))w(i)-bA(w^(i)ai(i) — 
1). By substituting C into the SG equation (30) and letting 
w^(i + l)ai(i -b 1) = 1, A is obtained as 

^ _ t/(i)x"(t)af (i) - (T?(t)af (i)ai(i) 

- imw) ■ ‘ ’ 

By substituting (31) back into (30) again, the weight update equation 
for LOCSME-SG is obtained as 


w(i-bl) = (I-/idi(i)ai(i)af (i))w(i) -y((I+ ) 

af (i)ai(t) 

t/*(I)x(i) - d-i(i)ai(i)). (32) 









The adaptive SG recursion circumvents a matrix inversion when 
computing the weights using (21), which is unavoidable in LOC- 
SME. Therefore, the computational complexity for computing 
the weights is reduced from 0{M^) in LOCSME to 0{M^) in 
LOCSME-SG (ISM^ + 25M). The proposed LOCSME-SG algo¬ 
rithm can be reproduced by using equations (9),(11),(12),(17),(22)- 
(28) and (32). Moreover, compared to existing RAB algorithms 
0 El [Hill [To] [n] and LOCSME which have a complexity equal or 
higher than 0{M^), LOCSME-SG has a greatly reduced cost. Com¬ 
pared with the approach in El which implements a low-complexity 
worst-case adaptive algorithm with a computational complexity of 
2M^ -f 7M, LOCSME-SG achieves a much better performance. 

5. SIMULATIONS 

A uniform linear array (ULA) of M = 12 omnidirectional sensors 
with a spacing of half wavelength is considered in the simulations. 
The desired signal is assumed to arrive at 0i = 10° while there 
are other two interferers impinging on the antenna array from di¬ 
rections 02 = 50° and 63 = 90°. The signal-to-interference ratio 
(SIR) is fixed at 20dB. We employ 100 repetitions to obtain each 
point of the curves and allow only one iteration performed per snap¬ 
shot. In our algorithm, the angular sector is chosen as [9i — 5°, 0i -I- 
5°] and the number of eigenvectors of the subspace projection ma¬ 
trix p is selected manually with the help of simulations. The pro¬ 
posed LOCSME and LOCSME-SG algorithms are compared with 
previously developed low-complexity standard SG algorithm and the 
RAB method which is based on the worst-case optimization Ql- 
We consider both coherent local scattering and incoherent local scat¬ 
tering scenarios for the mismatch and look at the beamformer output 
SINR in terms of snapshots with a maximum of i = 500 snapshots 
observed, or a variation of input SNR (-lOdB to 30dB) as shown in 
Figs. 1 and 2, respectively. 

5.1. Mismatch due to Coherent Local Scattering 

The steering vector of the desired signal affected by a local scattering 
effect is modeled as 


4 

a = p + (33) 

k=i 

where p corresponds to the direct path while h{9k){k = 1, 2,3,4) 
corresponds to the scattered paths. The angles 9k{k = 1,2,3,4) 
are randomly and independently drawn in each simulation run from 
a uniform generator with mean 10° and standard deviation 2°. The 
angles pkik = 1,2, 3,4) are independently and uniformly taken 
from the interval [0, 27r] in each simulation run. Notice that 9k and 
Pk change from trials while remaining constant over snapshots 0. 
We select p = 0.2, = 1, = 0.001, \q — 0.99, Ro = 101. 

The SINR performance versus snapshots and SNR of all the tested 
algorithms affected by coherent scattering is illustrated as in Figs. 1 
(a) and 2 (a). LOCSME-SG outperforms the other algorithms and is 
very close to the standard LOCSME. 

5.2. Mismatch due to Incoherent Local Scattering 

In the incoherent local scattering case, the desired signal has a time- 
varying signature and the steering vector is modeled by 

4 

a(i) ^ so{i)p + ^ Sk{i)h{0k), (34) 

k=i 


where Sk{i){k — 0,1,2, 3,4) are i.i.d zero mean complex Gaus¬ 
sian random variables independently drawn from a random gener¬ 
ator. The angles 9k{k = 0,1, 2, 3, 4) are drawn independently in 
each simulation run from a uniform generator with mean 10° and 
standard deviation 2°. This time, Sk{i) changes both from run to 
run and from snapshot to snapshot. We select /r = 0.1, pe = 5, 
(Te = 0.001, \q = 0.99, Ro = 501. The SINR performance versus 
snapshots and SNR is depicted in Figs. 1 (b) and 2 (b). Different 
from the coherent scattering results, all the algorithms have a certain 
level of performance degradation due to the effect of incoherent local 
scattering. However, over a wide range of input SNR, LOCSME-SG 
is able to outperform the standard SG and the low-complexity worst- 
case beamformers. 
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Fig. 1. SINR versus snapshots. 
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6. CONCLUSION 

The proposed LOCSME and LOCSME-SG algorithms only require 
prior knowledge of the angular sector of the desired signal and have a 
low complexity feature compared to prior art. Compared to the stan¬ 
dard SG beamformer and the low-complexity worst-case approach, 
LOCSME and LOCSME-SG have an outstanding performance in 
both coherent local scattering and incoherent local scattering cases. 
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